Radiography of the Earth's Core and Mantle with Atmospheric Neutrinos 



OO 
O 

o 

(N 
4=> 

0\ 



43 

Oh 1 
i ' 

CD 

43 



(N 
> 

in 
o 



o 



X 



M.C. Gonzalez-Garcia, 1,2 Francis Halzen, 3 Michele Maltoni, 4 and Hiroyuki K.M. Tanaka 5,6 

1 Institucio Catalana de Recerca i Estudis Avancats (ICREA), 
Departament d'Estructura i Constituents de la Materia, 647 Diagonal, E-08028 Barcelona, Spain 
''C.N. Yang Institute for Theoretical Physics, SUNY at Stony Brook, Stony Brook, NY 11794-3840, USA 
J Department of Physics, University of Wisconsin, Madison, WI 53706, USA 
■* Departamento de Fisica Teorica & Instituto de Fisica Teorica UAM/CSIC, 
Facultad de Ciencias C-XI, Universidad Autonoma de Madrid, Cantoblanco, E-28049 Madrid, Spain 
^Earthquake Research Institute, University of Tokyo, 113-0032 Tokyo, Japan 
6 Atomic Physics Laboratory, RIKEN, 351-0198 Saitama, Japan 

A measurement of the absorption of neutrinos with energies in excess of 10 TeV when traversing 
the Earth is capable of revealing its density distribution. Unfortunately, the existence of beams 
with sufficient luminosity for the task has been ruled out by the AMANDA South Pole neutrino 
telescope. In this letter we point out that, with the advent of second-generation kilometer-scale 
neutrino detectors, the idea of studying the internal structure of the Earth may be revived using 
atmospheric neutrinos instead. 
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The density profile and the shape of the Earth's core 
and mantle and their boundary (CMB) determine its geo- 
dynamo as well as the feeding mechanism of hotspots at 
the surface [jij . Knowledge of the CMB is derived from 
body-wave and free oscillation studies. The information, 
while more precise than what we can realistically expect 
from neutrino radiography in the near future, cannot re- 
duce ambiguities in our present model of the CMB asso- 
ciated with the fact that arrays of seismometers only pro- 
vide regional information, and that free-oscillation data 
only reveal one dimensional structure. The trade off 
among density, temperature, and chemical structure for 
body wave studies increases the uncertainty of the value 
for the density. For these reasons aspects of the global 
structure of the CMB region require confirmation. The 
study presented in this paper indicates that present neu- 
trino detectors have to be operated for 10 years to locate 
the CMB. IceCube will establish the averaged core and 
mantle density as a function of longitude thus providing 
the first independent global survey of the CMB region. 
We anticipate however that more precise global infor- 
mation on the CMB region will be obtained by longer 
observation periods or by future large scale neutrino de- 
tectors. Early studies of the possibility of doing neutrino 
tomography date back more than 25 years |2j. These 
proposed studying the passage of cosmic beams of high 
energy (HE) neutrinos through the Earth to diagnose its 
density. ^Alternatively, it was suggested to use accelerator 
beams [3] and to study the propagation effects through 
matter of oscillating neutrinos; for a review see Ref. [J]. 

The idea of neutrino tomography is straightforward: 
the Earth becomes opaque to neutrinos whose energy ex- 
ceeds ~ 10 TeV. The diameter of the Earth represents 
one absorption length for a neutrino with an energy ~ 
25 TeV. Such neutrinos are produced in collisions of cos- 
mic rays with nuclei in the Earth's atmosphere but, be- 



cause of the steeply falling energy spectrum of ~ E~ 3 7 of 
the atmospheric neutrino flux, such events are rare. The 
hope was that beams of cosmic neutrinos, likely to be as- 
sociated with the sources of the cosmic rays which reach 
energies of 10 8 TeV, would provide a plentiful source of 
neutrinos in the appropriate energy range. These would 
be detected by HE neutrino telescopes under develop- 
ment at the time. The cosmic beams would perform ra- 
diography of the Earth's interior as it moves relative to 
the comic source. In light of the recent development of 
successful and affordable technologies to build very large 
neutrino telescopes, we revisit the proposal ||. 

Neutrino telescopes detect the Cherenkov radiation 
from secondary particles produced in the interactions of 
HE neutrinos in deep water or ice. At the higher energies 
the neutrino cross section grows and secondary muons 
travel up to tens of kilometers to reach the detector from 
interactions outside the instrumented volume [6]. The 
construction of kilometer-scale instruments such as Ice- 
Cube at the South Pole and the future KM3NeT detec- 
tor in the Mediterranean, have been made possible by 
development efforts that resulted in the commissioning 
of prototypes that are two orders of magnitude smaller, 
AMANDA and ANTARES Q. Their successful tech- 
nologies have, in turn, relied on pioneering efforts by the 
DUMAND || and Baikal 0, as well as the Macro and 
Super-Kamiokande collaborations [13]. IceCube [lH is 
under construction and taking data with a partial ar- 
ray of 1320 ten inch photomultipliers positioned between 
1500 and 2500 meter, deployed as beads on 22 strings 
below the geographic South Pole. Its effective area al- 
ready exceeds that of its predecessor AMANDA by ~ 
one order of magnitude. The detector will grow by other 
14-18 strings in 2007-08 to be completed in 2011 with 
80 strings. 

AMANDA has observed neutrinos with energies as 
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high as ~100 TeV, at a rate consistent with the flux of at- 
mospheric neutrinos (ATM-z/s) extrapolated from lower 
energy measurements. It thus establishes limits on any 
additional flux of cosmic neutrinos in the energy range 
of interest for Earth tomography. These now reach be- 
low E^dN/ dE v < 1CT 11 TeVcm" 2 s" 1 sr" 1 for a diffuse 
flux [jj] and E 2 dN/dE v < 1(T 10 TeV cm" 2 s" 1 sr" 1 
From these results Standard Model physics is sufficient 
to establish that the event rates from cosmic beams per 
year in a future kilometer-scale detector are limited to 
~ 10 events from any particular source in the sky and 
less than ~ 100 from the aggregate of sources. Needless 
to say that the statistics is already uncomfortably small 
for a beam to be exploited for Earth tomography. 

Our main observation is that, with the growth of the 
detectors, the opportunity arises to exploit the ATM-y's 
that represent the background in the search for cosmic 
sources, as a beam for studying the Earth. The key point 
is that, the statistics for ATM-z/s in the 10 to 100 TeV 
energy range, is superior to those expected from any cos- 
mic sources detected in the future within the upper limits 
already established by AMANDA observations. Viewed 
from the South Pole, a uniform flux of ATM-z/s reaches 
the detector from the northern half of the sky; it will 
be modified in the 10 TeV energy region by its passage 
through the Earth. For instance, neutrinos from vertical 
to ~ 30 degrees have penetrated the core of the Earth 
before detection, whereas the ones detected at larger an- 
gles have traversed the mantle only. Establishing direct 
evidence for the transition from mantle to core will here 
be used as a benchmark to evaluate the technique. 

We will conclude that IceCube can directly observe 
the core-mantle transition at the 5cr level in 10 years. 
This evaluation is based on modeling of the HE ATM- 
i/'s that is, at present, still subject to uncertainties. We 
however establish that, under conservative assumptions, 
the transition can be observed at the 3<r level or above. 

We use the semianal ytic al calculation of IceCube event 
rates described in Ref. [14j . In brief, the expected number 
of z^-induced events (events arising from interactions 
can be evaluated similarly) in an exposure time T is: 
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is the differential flux in 



where d 2 <$> Vfi / {dE v d cos 
the vicinity of the detector after propagation in the 
Earth; more below. We use as input the neutrino fluxes 
from Honda [IH extrapolated to match the fluxes from 
Volkova [l6[ at higher energies. At the relevant energies, 
prompt v^s from charm decay are important. We intro- 
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FIG. 1: (a) Expected zenith angle distribution of ATM 
induced events in IceCube for different energy thresholds 
^.fin.min for the p R EM. 9 is the neutrino angle (which at 
these energies is collinear with the detected muon) as mea- 
sured from the vertical direction (upgoing-^ corresponding to 
9 = 180). (b) Ratio of the zenith angle distribution of ATM 
induced events in IceCube for different energy thresholds 
^fln,mm Qver correS p 0nc iing one for L2 cuts only. The 
shadow areas cover the angular size of the Earth core. 



duce them according to the recombination quark parton 
model [l7[ but consider alternative estimates also. 

d(TQ C / 'dE^EviE®) is the differential interaction cross 
section producing a muon of energy E®, After produc- 
tion the muon ranges out in the rock and in the ice 
surrounding the detector and loses energy to ionization, 
bremsstrahlung, e + e~ pair production and nuclear inter- 
actions. This is encoded in F(E®, E^ n 7 1) [l8| which rep- 
resents the probability that a muon produced with energy 
E® reaches the detector with energy E^ n after traveling 
a distance I. nx is the number density of nucleons in the 
matter surrounding the detector. 

The details of the detector are encoded in the effective 
area A® s for which we use the parametrization in Ref. 14 1 
describing the response of IceCube after the background 
rejection quality cuts referred to as "Level-2" (L2) cuts 
in Ref. [ll| . / m ; n = 300 m is the minimum muon track 
length required for the event to be detected. Effectively 
A° s vanishes for E% n < 100 GeV. 

In order to obtain d 2 (j) v /(dE v dco$9), one must ac- 
count for the simultaneous effects of oscillations and in- 
elastic interactions with the Earth matter which lead 
to the attenuation of the neutrino flux which are dif- 
ferent for v T 's and v^s. This can be achieved 
by solving a set of coupled evolution equations for the 
neutrino flux density matrix and for the muon and tau 
fluxes In practice, due to the steepness of the ATM- 



v spectra and the small value of the relevant Am 2 tm / E Vl 
d 2 (j) Vfi j \dE v dcos9) can be obtained from: 
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P^{E V ,L) exp{-X(^)[a N c(^) + *% G (E V )}}, (2) 

where P tlfJi (E v ,L = 2i?|cos#|) is the oscillation proba- 
bility. For E v > 1 TeV, P^ ~ 1. X(0) is the column 
density of the Earth, and R its radius: 

X{0)=N A / P e(VR 2 + z 2 + 2Rz cos 6) dz . (3) 
Jo 

iVyi is the Avogadro number, and p^(r) is the Earth mat- 
ter density assumed to be spherically symmetric. 

Equation ([2]) embodies the physics that makes Earth 
tomography with HE neutrinos possible. At sufficiently 
high energies, E v > 10 TeV, the attenuation fac- 
tor exp{— X(9)[a^c{E v ) + a cc (E u )]\ becomes relevant. 
Thus measuring Nev one can get information on p^,{r). 

In the left panel of Fig. [1] we show the expected zenith 
angle distribution of atmospheric ^-induced events in 
IceCube for different £* n > mm energy threshold as ob- 
tained using the Earth matter density profile of the Pre- 
liminary Reference Earth Model (PREM) In the 
PREM the Earth consists of a mantle extending to ra- 
dial distance r ~ 3000 km below the Earth surface and 
a core under it with a sharp core-mantle transition in 
density of about a factor 2. Thus neutrinos arriving with 
9 > 147 degrees (cos 8 < —0.84) will cross the core in 
their way to the detector. 

In the figure one notices, at sufficiently high energies, a 
reduction of the number of events for trajectories which 
cross the core resulting in a "kink" in the angular dis- 
tribution around 9 > 147. This feature is more clearly 
illustrated in the right panel of Fig. [TJ where we plot the 
ratio of the zenith angle distribution of events with ener- 
gies above E^ n - min divided by the number of events with 
no additional energy cut, which effectively corresponds 
to events with a threshold energy E^ n ' L2 ~ 100 GeV. 

Figure [T] illustrates the potential of doing Earth to- 
mography with the IceCube ATM-i/ samples. However 
one must realize that the angular dependence in the 
ratio shown in the right panel of Fig. Q] is not only 
due to Earth attenuation factor: there is an additional, 
Earth-independent, contribution from the variation of 
the zenith angle distribution of the fluxes with E v which 
does not cancel out in the ratio of events at different ener- 
gies. In principle this effect could be removed by compar- 
ing the ratio of upgoing (9 > 90 degrees) and downgoing 
events (9 < 90 degrees). In practice, the overwhelming 
atmospheric muon background makes the measurement 
of downgoing events impossible at these energies. 

In order to quantify the sensitivity of IceCube to the 
Earth density profile we study the ratio of observed 
events above a given energy threshold to the one expected 
for an Earth of equal mass as ours but with an homoge- 
neous matter distribution, phom = 3AiEarth/(47ri? 3 ): 

R = N^Ef->Ef^, cose, PPREM ) 
N^(E^>E h ^ mh \ cos 9, p hom ) 
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FIG. 2: Ratio of zenith angle distribution of expected events 
for the PREM over the expectations with an homogeneous 
Earth matter distribution for different values of the energy 
threshold of the events. The error bars in the figure show the 
expected statistical error in 10 years of IceCube. 
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TABLE I: Number of expected atmospheric ^-induced muon 
events in 10 years of IceCube operation in the different angu- 
lar bins and energy thresholds for the PREM. 



In Fig. [2] we show this ratio obtained by integrating the 
events in the numerator and denominator in 6 angular 
bins in cos (9, and for three values of the threshold en- 
ergy: E^ h2 ~ 100 GeV, 10 TeV, and 32 TeV. In this 
plot, trajectories crossing the core are contained in the 
most vertical bin. In the figure we also show the ex- 
pected statistical uncertainty <7 sta t,i> computed from the 
expected number of events in each angular bin in the 
PREM in 10 years of IceCube (see Table U). 

As expected, events with low energy threshold have 
no sensitivity to the Earth density and consequently the 
ratio for E^ n ' h2 ~ 100 GeV is practically constant and 
equal to 1. As E^ n ' mm increases the ratio becomes in- 
creasingly different from 1, reflecting the fact that the 
effect of the Earth matter profile becomes more evi- 
dent. The strategy is then obvious. One uses the mea- 
sured zenith angular distribution of the L2 event sam- 
ple as normalization to obtain the expectations for a 
constant density Earth at higher energies, N fJi (E^ n > 
£* n ' mm , cosf5, phom). By comparing the expectations 
with observation, one can quantify the sensitivity to the 
Earth matter profile. 

After normalizing to the observed L2 distribution 
residual theoretical uncertainties remain associated with 
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the predicted zenith angle distribution. They include sys- 
tematic effects in calibration, theoretical errors in the 
energy-angle dependence in the atmospheric fluxes due 
to the uncertainties in the K/tt ratio as well as in the 
contribution from charm (which are expected to be the 
largest at the relevant energies) , and the uncertainties in 
the neutrino interaction cross sections. Over the limited 
energy range relevant here, we conservatively account for 
those by introducing three systematic errors in the anal- 
ysis: an overall normalization error of 20%, and angular 
tilt uncertainty of 5% between horizontal and vertical 



events and an additional a s 



1% uncertainty due to 



uncorrelated systematics for each angular bin. With this 
we construct two simple \ 2 functions as 



Xnh = MiD 
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where Rf = (1 + 0.2£ norm ) (1 + 0.05 (cosfyfo-tut) and 
we have defined R c = R± and 
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with cr s 2 tatim 



X^taM and a s ys ,m = 4 x (!%) 2 - In 

»=2 

Eq. © Rf and R^ are the corresponding theoretical 
predictions including the normalization and tilt factors 
which minimize x na - ln ^ ne X 2 functions we have, con- 
servatively, not included the events in the most horizontal 
bin cos 9 > —0.17 where larger backgrounds from possi- 
ble remaining misreconstructed downgoing muons may 
be expected. In choosing the optimum energy threshold 
for this comparison, one has to take into account that, as 
the energy increases, the Earth matter profile becomes 
more evident but the statistics decreases and so does 
the achievable precision. For these simple observables, a 
compromise sensitivity is achieved for E^ n ' mm = 10 TeV. 

X nh quantifies the rejection power against the hypoth- 
esis of a homogeneous Earth density, while x 2 m gives the 
sensitivity to the specific difference in density between 
the core and the mantle. We find that, if no deviation 
from the PREM predictions is observed, in 10 years Ice- 
Cube can reject the homogeneity of the Earth with a 



the theoretical and systematic uncertainties are reduced 
to be below the statistical errors. Correspondingly, the 
difference in density between the core and mantle can be 
established with xl m = 9-22 (3cr-4.7cr). 

In summary, we conclude that IceCube will be able 
to measure the averaged core and mantle density with a 
significance that reaches, conservatively, 3<r in a decade. 
Our best guess is that a result can be obtained at the 5a 
level. 
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